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1. Introduction 

In these proceedings we summarize our results on the anomalous magnetic moment of the 
muon presented in [1] and [2]. This quantity provides a precision test of the Standard Model 
(SM), the electromagnetic (EM) [3], the weak [4] and the strong force. An overview of the dif- 
ferent contributions is given in Table 1 . Clearly the QED contributions are dominant, followed by 
the ones from QCD. It also evident that the strong contributions dominate the uncertainty of the 
SM prediction. Most relevant in that respect is the leading order hadronic vacuum polarisation 
(LOHVP) [5]. There is some hope that this uncertainty may decrease due to improved experi- 
mental input [6] as well as the efforts of the lattice community [7, 8, 9]. In the long run the most 
problematic contribution may therefore be hadronic light-by-light scattering (LBL), which cannot 
thus far be determined in a model-independent way. It has been considered in many approaches 
such as the Extended Nambu-Jona-Lasinio (ENJL) model [10], the Hidden Local Symmetry (HLS) 
model [11], vector meson dominance (VMD) approaches [12, 13], the non-local chiral quark model 
(NL;^QM) [14, 15], the chiral constituent quark model (;^CQM) [16], in holographic models [17] 
and Dyson-Schwinger Equations (DSEs) [18, 19, 1]. The lattice calculations of LBL are still at 
an exploratory stage [20]. The combined theory result quoted in Table 1 is a number that different 
groups agreed upon [21] by combining several hadronic models. For comparison, we list here the 
most recent experimental result from Refs. [22, 23] which reveals a discrepancy of about three 
standai^d deviations [5]; very recent results even indicate a significance ranging between 4.7 and 
4.9 a [24]. 

These results are exciting since they may provide a signal of new physics beyond the Standard 
Model. However, as we will argue in detail in the course of this work they should be taken with 
some caution. The present estimate of the hadronic light-by-light contribution including its poten- 
tial error has been derived from model calculations. Although much care has been invested in the 
error estimate, we believe that the inherent limitations of the models employed may have led to 
an overly optimistic value. Here we summarise our comparison between our approach and that of 
models for one of the more contentious contributions to LBL [1], in order to highlight this point. 
To this end we discuss the hadronic photon four-point function and summarize the DSE approach 
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QED 


116584718.1 ( 0.2) 


99.99390% 


00.00098% 


weak 


153.2( 1.8) 


00.00013% 


00.07910% 


QCD LOHVP 


6949.1(58.2) 


00.00596% 


82.69628% 


QCD HOHVP 


-98.4( 1.0) 


00.00008% 


00.02441% 


QCD LBL 


105 (26 ) 


00.00009% 


16.50391% 


Standard Model 


116591827.0(64 ) 


100% 


100% 


Experiment 


116592089 (63 ) 






Exp-Theo 


262 (89 ) 







Table 1: Standard Model contributions to the muon g — 2. 
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(a) (b) 



Figure 1: The two classifications of corrections to the photon-muon vertex function: (a) hadronic vac- 
uum polarization contribution to a^. The vertex is dressed by the vacuum polarization tensor Ylfiv', (b) the 
hadronic light-by-light scattering contribution to a^. 



used in this work in Sections 2 and 3. Then we shortly discuss the results for the leading order 
hadronic contribution in Section 4. In Section 5 we compare our method to the extended Nambu- 
Jona-Lasinio model. We explain similarities and differences and how these affect the results for 
the quark-loop light-by-hght contribution in Section 6. We conclude in Section 7. 



2. Basics 



On mass-shell, the muon-photon vertex can be written in terms of two momentum dependent 
form-factors 




u{p') 



2m,, 



u{p). 



(2.1) 



where p and p' are the muon momenta, q is the photon momentum and (7„jg = 2[7ai7/3]- The 
anomalous magnetic moment is defined in the limit of vanishing photon momentum, q^, as 



g-2 



F2i0). 



(2.2) 



The HVP contribution, shown in Fig. 1(a) is defined within the DSE approach as the hadronic 
photon self energy 



n^v = 'vv^^^^^Vv, 



(2.3) 



which is itself defined in terms of the full quark and the quark-photon vertex. This representation 
is exact, up to truncations of vertex and propagator. This quantity is well described within our 
approach as is detailed in section 4. 



The LBL contribution to this vertex is depicted in Fig. 1(b). This hadronic photon four-point 
function can be split into several parts, which we organize by the number of quark T-Matrices 
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Figure 2: The quark loop contribution to the muon g — 2. The quarks and vertices are fully dressed, 
involved 




(2.4) 



The pre factors indicate the number of permutations that must be computed. The second term, 
containing just one T-matrix, has been considered in Refs. [18, 19] by performing a pseudoscalar 
resonant expansion. We found good agreement with corresponding results from hadronic models 
[10, 25, 26, 27, 28, 12, 13, 29, 21, 14, 30, 16]. The ellipsis contains, for example, contributions that 
can be written as charged pion loops including also the pion polarizabilities. While in earlier works 
these contributions have been argued to be small, a recent work suggests substantial contributions 
to LBL which need to be evaluated carefully in the future [31]. Note, however, that within the ap- 
proximation of QCD, we are working with (see Sec. 3) the representation in Eq. (2.4) is consistent 
with the exact representation Eq. (2.3), a fact that is interesting when it comes to considerations 
about systematic uncertainties. 

Later we concentrate on the first part of Eq. (2.4) referred to here as the quark-loop. It consists 
of four fully dressed quark propagators as well as four fully dressed quark-photon vertices. Its 
contribution to is shown explicitly in Fig. 2. The focus of these proceedings will be in the 
transverse parts of these fully-dressed vertices, which in principle contain contributions from vector 
mesons. It has been shown that these contributions generate a strong suppression of the quark-loop 
contribution to in the ENJL model [10]. Below we will argue that this suppression is in fact a 
model artefact attributable to the contact interaction featured therein. 

3. Method 

First, let us briefly review the needed Dyson-Schwinger equations and their truncation. The 
first ingredient for the calculation of the quark-loop is the quark propagator, S{p), whose inverse is 
described by the DSE 

S-\p) =Z2So' +g^ZifydkrS{k)r''{k,p)D^y{q), (3.1) 

To solve this equation we require the dressed gluon propagator D^y{p) and the dressed quark- 
gluon vertex Yv{k,p). In the rainbow truncation employed here we use a combined ansatz for both 
of these quantities retaining only the 7v component of the vertex. This truncation has been proposed 
in Refs. [32, 33] and has achieved phenomenological success in describing meson properties such 
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as masses, decay constants and electromagnetic form factors [32, 33, 34, 35, 36] as well as various 
baryon properties [37, 38, 39, 40]. 

The second ingredient for the calculation of the quark-loop is quark-photon vertex, F^, deter- 
mined from the inhomogeneous equation 

[r^(P,fc)]„ = Zi7^ - I dq[S{q+)T^{P,q)S{q-)]utKtuAk-q), (3-2) 

where Ktu.rs {k — q) is the quark-antiquark interaction-kernel. To satisfy the vector- and axial-vector 
Ward-Takahashi identities we choose the ladder-exchange kernel, consistent with rainbow trun- 
cation of the quark DSE above. An important consequence of this truncation is the dynamical 
generation of vector meson poles in the transverse components. Thus the tenets of vector me- 
son dominance models (VMD) are automatically encoded on the level of interacting quarks and 
gluons. Note furthermore that the vertex from Eq. (3.2) and the quark from Eq. (3.1) fulfill the 
Ward-Takahashi identity (WTI) 

iP^T^,{P,k) = S-\k^)-S-\k+). (3.3) 

Both properties are crucial for the realistic description of electromagnetic properties, e.g. the pion 
charge radius [34] . Their relevance to the leading order hadronic vacuum polarisation contribution 
to has been highlighted in Ref. [2], see Sec. 4. We will see below also the importance of the 
correct description of the transverse structures as regards the hadronic LBL contribution. From this 
we will argue that the twenty percent error quoted in Table 1 is too optimistic since it is obtained 
from approximations that implicitly miss information that we later find to be important. 



4. The leading order QCD contribution 

In the following we summarize our results for the leading order hadronic contribution, the 
leading order hadronic vacuum polarisation published in [2]. This quantity is a very important test 
for any approach that aims at a quantitative prediction of the much more complicated LBL. In short, 
using the quark [Eq. (3.1)] and the quark-photon vertex [Eq. (3.2)] we calculate this contribution 
via Eq. (2.3). Our result is 

4^^™^ = (6760-7440) xlO-i\ 

where the given range reflects our systematic model uncertainty, see [2] for details. Our result is 
close to the rather model independent analysis from Ref. [5] , i.e. 6949 x 10^^'. This is a deviation 
of less then ten percent. Although this precision is not a challenge for the model-independent HVP 
predictions, a con^esponding precision for LBL would be a great breakthrough. Furthermore, in 
Ref. [2] we presented a determination of the Adler function, that can be obtained from the hadronic 
vacuum polarisation tensor, c.f. Eq. (2.3), that shows good agreement with model-independent re- 
sults from dispersion relations on all scales. What is especially interesting in the context of these 
proceedings is that we found the transverse part of the quark-photon vertex to yield the dominant 
contribution of ~ 80% to a^^^^ . Since this part contains a dynamically generated vector meson 
bound state, this explains why VMD approximations work reasonably well for HVP, see e.g. [41]. 
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Figure 3: The LOVHP contribution to as a function of the p mass. We compare our results for the 
Nf = 2 (black) and Nf = 2 + 1 (red) flavour cases to lattice results. The grey bands represent our numerical 
uncertainty. 

Models that miss these physics have to use a constituent quark mass much lower than expected 
[16, 42]. 

In Ref. [2] we also compared our results against lattice QCD. Here, we present an updated 
comparison including more recent lattice data in Fig. 3. Both, the lattice and the DSE calculations 
are performed at various current quark masses. In order to provide a gauge and scheme invariant 
comparison between the approaches we have plotted the results against the corresponding varying 
mass of the p-meson. In the figure, we compare against results from Refs. [7, 8, 43] for the cases 
of Nf = 2 (black) and Nf = 2 + I (red) flavours. We find agreement within error bars, showing 
nicely that the DSEs have the correct physics included and that the p meson mass is the important 
scale here. 

It is interesting, that while we confirm the simple picture of VMD for the case of the leading 
hadronic contribution, the picture will be very different in the case of the LBL contribution. This 
is the subject in the following. 

5. Comparison: DSE vs. ENJL 

To explicate the differences between our approach and others, we make a comparison between 
the DSE truncation used in this work and the ENJL model. The first striking difference lies in 



6 



Role ofVMD in hadronic LBL contribution to the muon g — 2 



Christian S. Fischer 



the fact that DSEs have a smooth ultraviolet limit, wherein at large momenta we connect with 
perturbation theory. ENJL on the other hand is non-renormalisable and features an effective cut- 
off scale on the order of a GeV. This necessitates the splitting of the photon four-point function 
into a high- and low-energy region that is not unique due to their being two independent scales. 
Additionally, the dressings of the ENJL model feature trivial momentum dependencies due to its 
contact interaction. In this regard the DSE approach is quite different. In what follows we will 
explore this is some detail; further information can be found in Ref. [1]. 

Primarily we will focus on the two flavour case as this is sufficient to make our point. However, 
we will quote Nf = 4 results for completeness in the later sections. 

5.1 The ENJL perspective 

The discussion here is based on Refs. [10, 44], where further details may be found. In the ENJL 
model, the quark propagator assumes the form of a momentum independent constituent quark 

^ENJL(/') = -'y+^> (5-1) 

with a typical mass of M 300 MeV. The wave function renormalisation, Zy , is unity. The quark- 
photon vertex is given by the bubble sum 



(5.2) 



which can be resummed as a geometric series ^„ Bubble" = 1/(1 — Bubble) . This simple behaviour 
is due to the effective contact interaction that decouples the loop integrals. Consequently, the 
result depends only on the total photon momentum Q and has no dependence on the relative quark 
momentum. Explicitly the vertex has the form 

pENJL_ T 

which contains the bare vertex 7^ and the leading transverse structure = {5^v — Q^Qv/Q^)Yv 
The dressing of this transverse part is given here in the VMD limit of the ENJL model [44] where 
for two flavours My is identified with the p-mass. Using the transversality of the hadronic photon 
four-point function with respect to its photon legs, this vertex can be reduced to y^My /{Q^ +My) 
which highlights the strong suppression induced by the VMD contribution to the vertex. Note that 
this model is consistent with the Ward-Takahashi identity, with the 7^ component the constrained 
gauge-part and 7^ the transverse part. 

5.2 The DSE perspective 

Now we investigate the corresponding objects in the case of DSEs. The full inverse quark 
propagator is given by 

S^Up)=Z-f.\p^){-i/ + M{p^)), (5.4) 
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which contains two momentum dependent dressing functions Zf{p^) and M{p^), which correspond 
to the wavefunction renormalisation and the mass function, respectively. The quark-photon vertex 
is given as a sum of gluon ladders 



S 



S 



(5.5) 



which, as opposed to Eq. (5.3), can not be resummed as a trivial series and must instead be solved 
numerically. In general, this vertex can be decomposed into twelve tensor components 



(G, ^) = £ A (2, k)l^ + £ t(') (G, fc) T, 



^(0 



(5.6) 



(=1 



j=i 



where k is the relative quark momentum and Q the total photon momentum. Note that the first part, 
containing the A^'', is also called the Ball-Chiu vertex r^*- [45]. The vector meson bound state does 
also appear in the transverse vertex structure in the DSE/BSE case. This is most easily pictured in 
the form of a simple fit [34] to the numerical results of the quark-photon vertex 



-7/ 



r (o'Nv fv 



(5.7) 



which consists of the non-transverse Ball-Chiu part, T^, and the leading transverse structure coitc- 
sponding to = 7^. For the parameters we find reasonable agreement with the numerical solution 
with ft) = 0.66GeV, a = 0.15 and Nv/v/My = 0.152. Note that, as in the ENJL model, we have 
in Eq. (5.7) a part that is given via the WTI (the BC vertex) and a transverse part. We wish to 
emphasize once more, that the vector meson pole appearing in the transverse part is generated 
dynamically in the DSE and ENJL approaches. 

5.3 Differences between DSEs and ENJL model 

We already stated that one of the chief differences between the DSE and ENJL approach lies 
with the momentum dependencies of the propagators and vertices. We compare the quark prop- 
agators in Fig. 4, where one sees a non-trivial momentum dependence in both the wave function 
renormalisation and the mass function for the DSE. In particular, the mass function features the 
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Figure 4: Comparison of the mass function M and the wavefunction Zf for DSE and ENJL quarks. 
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expected resumed one-loop perturbative behaviour at large momenta and smoothly connects to the 
infrared domain where a constituent-like mass is dynamically generated. This rapid accumulation 
of mass is consistent with results of the operator product expansion [46]. Indeed, the same qual- 
itative behaviour is seen on the lattice, see Ref. [47] for a comparison. In contrast, the dressing 
functions in the ENJL model are just constants and hence miss essential features present in QCD. 

Next we compare the approximate functional forms of the quark-photon vertex in the ENJL 
model, Eq. (5.3) and DSE approach, Eq. (5.6). The dressing of the 7^ part is shown in Fig. 5. Con- 
sistent with the Ward-Takahashi identity it is unity in the ENJL case due to the quark wavefunction 
being trivial there. In contrast, the corresponding vertex dressing in the DSE approach behaves as 
A^') ~ l/Zf, owing again to Ward identities. A consequence of the wavefunction being screened 




Figure 5: The leading A^^' component of the quark-photon vertex constrained by the Ward-Takahashi iden- 
tity of Eq. (3.3). We show two slices with relative, k, and total, Q, momentum set to zero, respectively. The 
constant dressing corresponds to the ENJL model. 



for a large part of the momentum regime is that A^^^ > 1 leading to an enhancement. However, 
since each vertex can be paired with a quark propagator that features a screening ~ Zf these effects 
essentially cancel in the quark- loop contribution to . Thus, quantitative differences between the 
ENJL model and the DSE approach owing to non-trivial wavefunction renormalisation and the 
A*-^)/^ part of the vertex are expected to be small. This is confirmed in explicit calculations below. 

Now we consider the transverse part of the vertex, T^^' = 7j[. We compare the fit function 
of Eq. (5.7) with the full numerical solution, of the vertex BSE in Fig. 6 for a momentum 
slice in which the relative momentum between the two quarks is vanishing. We see that at small 
photon momentum, Q, qualitative agreement between the ENJL and DSE calculations. Since the 
ENJL model is not ultraviolet complete, it is not surprising that large differences arise at lai^ge 
values of the photon momentum. Fortunately, under the integral such a momentum region is not 
weighted strongly in the determination of and so this deviation is not important. This situation 
is drastically changed when one takes into account the relative momentum between the two quarks, 
i.e. 7^ 0. We see from the fit to the DSE calculation of the vertex, Eq. (5.7) that the transverse 
part of the vertex is highly suppressed for > co^. The ENJL model presents no such dependence. 
The quantitative impact of this damping will be discussed in the next section. 
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Figure 6: The dependence of the dominant transverse dressing t''^ on the space-like photon momentum 
is shown for the explicit solution of the quark-photon vertex BSE, Eq. (3.2), the fit to this given in Eq. (5.7), 
and the transverse part of the ENJL vertex, Eq. (5.3). 



6. The hadronic light-by-light contribution 



We are now ready to turn the qualitative discussion in the last section into quantitative state- 
ments. To this end we give a summary of the results of Ref. [1] on the contribution to Off^' quark-ioop 
using various quark- and vertex-dressings. For some quantities we calculate the averages of dress- 
ing functions as the arithmetic mean from samples of the VEGAS Monte-Carlo routine [48] we 
employ for a^^' ^"^'^"'°°p \Ye consider the case of Nf = 2 flavours only. The remaining flavours 
constitute a correction of the order of ten percent and don't change the general picture that emerges 
from this comparison. 

6.1 Influence of the quark and gauge part of the vertex 

Here, we explicitly show that the suspected cancellation between the wavefunction renormal- 
isation, Zf, of the quark propagator and the leading gauge part of the quark-photon vertex does, 
indeed, take place in the calculation of the quark- loop contribution to a^. Since the dressing of 
the quark photon vertex is directly related via the Ward-Takahashi identity to Zf it suffices to per- 
form the calculation in full, and with Zf = \jZ'i . We display the results in Table 2. As is evident 







M,Z/,A(') all dynamical 




100 


M dynamical, Zf = 1 /Z2, 


A(i)=Z2 


102 


M = 0.2GeV, Zf= I/Z2, 


A(i)=Z2 


104 



Table 2: The quark-loop contribution to hadronic LBL. Compared are our result for full quark 

propagator and gauge part of the vertex with an approximation where the momentum dependence of the 
quark wave function is neglected. In addition we show the result where also the quark mass is replaced by a 
constant M = 0.2 GeV. 

'The Z2 is a necessary renormalisation constant to ensure multiplicative renormalizability 
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by the similarity of the two results, this clearly supports our qualitative discussion regarding the 
cancellation. 

In the process of performing these integrals, we also determined the average mass probed to be 
« 0.2 GeV, as weighted by the Monte-Carlo integration. While this finding explains why 
constituent quark models such as Ref. [16] and [49] must necessarily use very light constituent 
quark masses in order to achieve reasonable results, it also shows that this definition of effective 
mass is process dependent, contrary to the running mass function defined through the quark DSE, 
which is universal. The projection onto a certain process is only introduced through the process 
dependent Feynman-diagrams which includes the quark and thus the mass function. 

We wish to emphasize that the X^^^Jii P^t oi^ly oris of three terms in the gauge part of 
the vertex. Although the corresponding result is already approximately gauge invariant [1], the 
functions A^^/^' from Eq. (5.6) have to be included for principal reasons. Unfortunately the cor- 
responding numerics turns out to be delicate [19], such that more work is needed to resolve this 
issue. 

6.2 Influence of the vertex: transverse part 

Now we investigate the impact of the leading transverse part of the vertex, 7^, on the quark- 
loop contribution to . 

Since the ENJL model does not feature a relative momentum dependence in its vertex dress- 
ings, we compare the ENJL result for the transverse vertex with our DSE determinations where we 
force the condition = 0. From the left-panel in Table 3 we see agreement between the results 
which highlights that the differences at large momenta in the fit functions are largely irrelevant. 

In the right-panel of the same table, we show the corresponding results for when the relative 
quark momentum is taken into account. An appropriate function is also introduced to the ENJL 
parametrisation 

with CO defined above. Again, we see a degree of parity between the models but note that there is 
more than a factor of two difference as compared to the results with restricted momentum depen- 
dence. 



Vertex Dressing, = 




Vertex Dressing, Z:^ 7^ 






43 


7^.^(1) + 7^Tenjl/(^') 


103 




43 


7.A(') + 7^^:^ 


105 


7.AW+7MVc 


41 




96 



Table 3: Leading gauge part and leading transverse vertex component, with dressing functions from the 
ENJL model, VMD like fit from DSE/BSE (t^J'), and from an explicit calculation of the quark-photon DSE 
'^^caic-'- Results are shown without (left table), and with (right table) the inclusion of a dependence on the 
relative momentum. 

Our most sophisticated vertex construction is the leading gauge part of the vertex plus the 
leading transverse part, corresponding to the last line of the right Table 3. Repeating this calculation 
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for the case of A^j = 4 flavours we obtain 

a;''^'^''^^=' = (107±2)xl0-", (6.2) 
where the error is purely statistical. 

7. Summary and conclusions 

We summarised the detailed comparison of the DSE framework to the ENJL model with re- 
spect to the quark loop part of the hadronic light-by-hght contribution to the anomalous magnetic 
moment of the muon, presented in Ref. [1]. We found similarities but also important dif- 

ferences. Our main focus was on the influence of the momentum dependence of dressing functions 
that are not present in the simplified treatment within the ENJL model. These are the quark mass 
and wave functions as well the leading non-transverse and transverse dressing functions of the 
quark-photon vertex. Since some of the momentum dependencies are related via a Ward-Takahashi 
identity they cancel approximately in the calculation of Off^''^'. Other important differences, how- 
ever, remain. One is the impact of the running quark mass function, which turns out to be equivalent 
to an effective constituent quark mass of M « 0.2 GeV, considerably smaller than typical values 
considered within the ENJL model. This should, however, not be understood as a general definition 
of a constituent mass, but rather a way how to make a complicated integration simpler to get some 
intuitive insight into the physics and numerics at work. The other important difference concerns 
the impact of the transverse parts of the vertex. In the ENJL model, the dependence of this part 
of the vertex on the relative momentum k between the quarks is neglected resulting in a large de- 
crease of Off^''^' due to transverse contributions. Within the DSE approach we have shown that this 
suppression is dramatically reduced when the full momentum structure of the vertex is taken into 
account. 

In Table 4 we compare our current results for the different contributions to a^^^ with the 
ENJL-results of Ref. [10], the current 'standard' result of Ref. [21] (summarising calculations from 
different sources) and the values from a recent estimate within a chiral quark model [16], which 
has been advocated as a qualitative reference. Compared to the ENJL-model we nicely agree in 



[10-11] 


LBLXko) 


LBL,m,n') 


LBLmI 


LBL.other 


„LBL 
"/i 


ENJL,Ref.[10] 


85 (13) 


21(3) 


-23 (16) 


83 (32) 


standard,Ref.[21] 


116(13) 


-11(13) 


105 (26) 


;t;QM,Ref.[16] 


68 




82 




150 


DSE, Ref. [1] 


58(1) 


23(1) 


107 (2) 




188 (4)(90) 



Table 4: Selected partial and total results for . In the first line we display the ENJL result, the second 
line displays the current 'standard result' (numerical and systematic error), the third line stems from a recent 
estimate within a chiral quark model [16] (evaluated at a constituent mass of M = 240 MeV), whereas the 
fourth line shows our results so far (numerical error only, with a guess of the systematic error of the total 
result). Cells with a dash mark contributions which have not (yet) been calculated in the coiTesponding 
approach. The column 'other' sums up contributions from scalar and axialvector meson exchange as well as 
contributions from dressed pion loops. 
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the meson exchange contributions a'jf^'^'^^ and a'jf^'^^'^ ^ but note a drastic enhancement of the 
quark-loop contribution. We have explained this difference above. The same arguments apply to 
the results of the summary of Ref. [21] in the second line. The comparison with the chiral quark 
model (;t;QM) in the third line is interesting. They have applied a constituent mass in the range 
of M = 230 — 250 MeV. As we have argued above, our results with full momentum dependence 
are reproduced by an approximation using an average quai^k mass function of (M) = 200 MeV, i.e. 
even lower. For such a small value the results of Ref. [16] become smaller in the pion exchange 
channel but larger in the quark-loop contribution resulting in a total of 168 x 10^^^ in very good 
agreement with our full results. 

Small constituent quark masses of order ~ 200 MeV have been found to be necessary in 
various approaches using simple models, see e.g. Refs. [50, 42, 49]. As argued above, our findings 
can explain this on a much deeper level of sophistication. It has to be mentioned, however, that 
not all models agree on whether the quark-loop- and meson-exchange-contribution are two ways 
to describe the same contribution, or whether these have to be added. From the viewpoint of QCD 
diagrams it seems clear that these are completely separate contributions, see e.g our Ref [19]. 

We conclude from our study that the standard value a'jf^ = 105(26) x [10^^'] used in current 
evaluations of the anomalous magnetic moment of the myon [5, 24] may be too small concerning 
its central value and is probably much too optimistic in its error estimate. We emphasize once 
again the importance of taking explicit momentum dependencies into account, a task that the DSE 
framework is well capable for. Our current results for a'jf^ are preliminary, since we neglected 
contributions from the transverse and non-transverse parts of the quai^k-photon vertex in the quark- 
loop. We are working on this issue. We are also working on an explicit calculation of contributions 
from the charged pion-loop which, according to the study of [31] may be more important than 
previously thought, see also the contribution of Ramsey -Musolf to these proceedings [51]. 

To compete with the expected enor bars of the forthcoming experiments we need to improve 
the theoretical error of a'jf^ into the twenty, or even better into the ten percent region. Within the 
DSE framework this level has been reached in the case of hadronic vacuum polarisation [2], where 
nice agreement with corresponding lattice calculations has been obtained. For hadronic light-by- 
light we propose a similar strategy. We believe that only the systematic comparison of different 
methods Uke effective models, the DSE framework and lattice gauge theory provides the potential 
for a reliable calculation of a^^^. 

Acknowledgments 

This work was supported by the DFG under contract FI 970/8-1 and the Austrian Science Fund 
FWF under project M1333-N16. 

References 

[1] T. Goecke, C. S. Fischer and R. WilUams, Phys. Rev. D in press, arXiv: 1210.1759 [hep-ph] . 

[2] T. Goecke, C. S. Fischer, and R. Williams Phys.Lett. B704 (201 1) 21 1-217, arXiv :1107.2588 
[hep-ph] . 

[3] T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio arXiv: 1205 . 5370 [hep-ph] . 

[4] A. Czarnecki, W. J. Marciano, and A. Vainshtein Phys.Rev. D67 (2003) 073006, 

arXiv:hep-ph/0212229 [hep-ph]. 



13 



Role ofVMD in hadronic LBL contribution to the muon g — 2 



Christian S. Fischer 



[5] K. Hagiwara, R. Liao, A. D. Martin, D. Nomura, and T. Teubner J.Phys.G G38 (201 1) 085003, 

arXiv: 1105 . 3149 [hep-ph]. 

[6] F. Jegerlehner and A. Nyffeler Phys.Rept. 477 (2009) 1-1 10, arXi v :0902.3360 [ hep-ph ] . 

[7] X. Feng, K. Jansen, M. Petschlies, and D. B. RenneT Pliys.Rev.Lett. 107 (2011) 081802, 

arXiv: 1103 . 4818 [hep-lat]. 

[8] P. Boyle, L. Del Debbio, E. Kerrane, and J. Zanotti arXiv: 1107 . 1497 [hep-lat] . 

[9] M. Delia Morte, B. Jager, A. Juttner, and H. Wittig JHEP 1203 (2012) 055, arXi v :1112.2894 

[hep-lat ] . 

[10] J. Bijnens, E. Pallante, and J. Prades Nucl.Phys. B474 (1996) 379-420, arXiv : hep-ph/ 95 1 1388 
[hep-ph] . 

[11] M. Hayakawa, T. Kinoshita, and A. Sanda Phys.Rev.Lett. 75 (1995) 790-793, 

arXiv:hep-ph/9503463 [hep-ph]. 

[12] M. Knecht and A. Nyffeler Phys.Rev. D65 (2002) 073034, arXiv : hep-ph/ 111058 
[hep-ph] . 

[13] K. IVIelnikov and A. Nzmshidn Phys.Rev. D70 (2004) 113006, arXiv: hep-ph/ 03 1222 6 
[hep-ph] . 

[14] A. E. Dorokhov and W. Broniowski Phys.Rev. D78 (2008) 07301 1, arXiv:0805.0760 
[hep-ph] . 

[15] A. Dorokhov, A. Radzhabov, and A. Zhevlakov arXiv: 1204.3729 [hep-ph] . 

[16] D. Greynat and E. de Rafael, ///fiP 1207 (2072) 020 arXiv: 12 04 . 302 9 [hep-ph]. 

[17] L. Cappiello, O. Cata, and G. Ambm&io Phys.Rev D83 (2011)093006, arXiv: 100 9 . 1161 
[hep-ph] . 

[18] C. S. Fischer, T. Goecke, and R. WilHams Eur.Phys.J. A47 (201 1) 28, arXiv :1009.5297 
[hep-ph] . 

[19] T. Goecke, C. S. Fischer, and R. Wilhams Phys.Rev D83 (2011) 094006, arXiv: 1012 . 3886 
[hep-ph] , Erratum-ibid. D86 (2012) 099901. 

[20] T. Blum, M. Hayakawa and T. Izubuchi, arXiv: 1301 . 2607 [hep-lat] author = "Blum, T. and 
Hayakawa, M. and Izubuchi, T.", title = "Hadronic corrections to the muon anomalous magnetic 
moment from lattice QCD", year = "2013", eprint = "1301.2607", archivePrefix = "arXiv", 
primaryClass = "hep-lat", SLACcitation = " 

[21] J. Prades, E. de Rafael, and A. Vainshtein arXiv: 0901 . 0306 [hep-ph] . 

[22] IVIuon G-2 Collaboration, G. Bennett et al. Phys.Rev D73 (206) 072003, 

arXiv:hep-ex/0602035 [hep-ex]. 

[23] B. Roberts Chin.Phys. C34 (2010) 741-744, arXiv: 1001 . 28 98 [hep-ex] . 

[24] IVI. Benayoun, P David, L. DelBuono and F. Jegerlehner, arXiv; 1210.7184 [hep-ph]. 

[25] J. Bijnens, E. Pallante, and J. Pmdes Nucl.Phys. B626 (2002)410^11, arXiv: hep-ph/ 01 122 55 
[hep-ph] . 

[26] M. Hayakawa, T. Kinoshita, and A. Sanda Phys.Rev D54 (1996) 3137-3153, 

arXiv:hep-ph/9601310 [hep-ph]. 



14 



Role ofVMD in hadronic LBL contribution to the muon g — 2 



Christian S. Fischer 



[27] M. Hayakawa and T. Kinoshita Phys.Rev. D57 (1998) 465^77, arXiv: hep-ph/9708227 
[hep-ph] . 

[28] M.HayakawaandT.KinoshitaarXiv:hep-ph/0112102 [hep-ph]. 
[29] A. Nyffelerfo5CD09 (2009) 080, arXiv: 0912. 1441 [hep-ph]. 

[30] A. Dorokhov, A. Radzhabov, and A. Zhevlakov Eur.Phys.J. C71 (2011) 1702, arXiv : 1103 . 2 042 
[hep-ph] . 

[31] K. T. Engel, H. H. Patel and M. J. Ramsey-Musolf, Phys.Rev. D86 (2012) 037502, arXiv: 1201.0809 
[hep-ph] 

[32] R IVIaris and C. D. Roberts Phys.Rev. C56 (1997) 3369-3383, arXiv : nucl-th/97 0802 9 
[nucl-th] . 

[33] R Maris and R C. Twdy Phys.Rev. C60 (1999) 055214, arXiv : nucl-th/ 9905 05 6 
[nucl-th] . 

[34] R IVIaris and R C. Twdy Phys.Rev C61 (2000)045202, arXiv:nucl-th/9910033 
[nucl-th] . 

[35] D. Jarecke, R Maris, andR C. Tw\dy Phys.Rev. C67 (2003) 035202, arXiv : nucl-th/02 08 01 9 
[nucl-th] . 

[36] R Maris andR C. Tandy Phys.Rev C65 (2002)045211, arXiv : nucl-th/02 1017 
[nucl-th] . 

[37] G. Eichmann, R. Alkofer, A. Krassnigg, and D. Nicmorus Phys.Rev.Lett. 104 (2010) 201601, 

arXiv: 0912 . 2246 [hep-ph]. 

[38] G. Eichmann Phys.Rev D84 (201 1) 014014, arXi v : 1 1 4 . 4 5 5 [ hep-ph ] . 

[39] G. Eichmann and C. S. Fischer, Ewr. P/j>x 7. A 48 (2072) 9, arXiv: 1111 . 2614 [hep-ph]. 

[40] H. Sanchis-Alepuz, G. Eichmann, S. Villalba-Chavez, and R. Alkofer P/zys./^ev. D84 (2011) 096003, 
arXiv :1109.0199 [hep-ph] ; H. Sanchis-Alepuz, R. Alkofer, G. Eichmann, and R. WilHams 
Po5 QCD-TNT-II (2011) 041, arXiv: 1112 . 3214 [hep-ph] ; H. Sanchis-Alepuz, R. Alkofer, 
andR. WilHams Po5QNP2012 (2012) 112, arXiv: 120 6 . 65 99 [hep-ph]. 

[41] J. S. Bell, E. de Rafael, Nucl. Rhys. Bll (1969)611-620. 

[42] R. Boughezal and K.MelnikovP/m.Leff. 6704(2011) 193-196, arXiv: 1104 . 4510 [hep-ph]. 

[43] M. Delia Morte, B. Jager, A. Juttner et al, [arXiv: 101 1.5793 [hep-lat]]. 

[44] J. Bijnens and J. Prades Z.Phys. C64 (1994)475-494, arXiv: hep-ph/ 94 032 33 [hep-ph] . 

[45] J. S. Ball and T.-W. Chiu Phys.Rev D22 (1980) 2542. 

[46] C. S. Fischer, J. Rhys. G 32 (2006) R253 [hep-ph/0605173]. 

[47] C. S. Fischer, D. Nickel, and J. Wambach Phys.Rev D76 (2007) 094009, arXiv:0705.4407 
[hep-ph] . 

[48] T. Hahn Comput.Phys.Commun. 168 (2005) 78-95, ai-Xiv:hep-ph/0404043 [hep-ph]. 

[49] R Masjuan and M. Vanderhaeghen, arXiv:1212.0357 [hep-ph]. 

[50] A. A. Pivovarov, Phys. Atom. Nucl. 66, 902 (2003) [Yad. Fiz. 66, 934 (2003)] [hep-ph/01 10248]. 

[51] M. J. Ramsey-Musolf, these proceedings. 



15 



